Pakistan Suicide Bombing Attacks (1995 - 2016) 


01. Uploading the Dataset:- 


import pandas as pd 


bombari df = pd.read csv("PakistanSuicideAttacks Ver 11 (30-November-2017).csv",encoding-'cp1252') 


bombari_df 
^ Blast " 
si Date same Day Holiday time City 
Date Type 
Type 
Sunday- pes 
0 1 November al-THaany Holiday Weekend NaN Islamabad 
19-1995 1416 AH 
Monday- 10 " 
1 2 November  SHa'baan Mie NaN NaN Karachi 
6-2000 — 1421 AH ay: 
Wednesday- 25 safar Working A " 
2 3 May 8-2002 — 1423 AH Day NaN 7:45 AM Karachi 
E 3 Raby’ al- - Wn 
2) Steet THaany vene NaN La Karachi 
j 1423 AH Y 
" 4 Jumaada ` 
4 5 ka al-awal working NaN NaN Quetta 
7 1424 AH Y 
Thursday- 13 " 
491 492 October-05- MuHarram ides NaN NaN Quetta 
2017 1439 AH ay 
Thursday- 27 , 
492 493 October-19- MuHarram MED NaN 8:00 AM Quetta 
2017 1439 AH gu 
Thursday- A 
493 494  November- A ieu NaN 7:00AM Quetta 
9-2017 0 ay 
Friday- 5 Rabi Al- ; 
494 495  November- Awwal HENDE NaN 7:00AM Peshawar 
24-2017 1439 AH ay 
Saturday- 6 Rabi Al- Worki 
495 496 November- Awwal is NaN 7:00AM Quetta 
25-2017 — 1439 AH ay 


496 rows x 26 columns 


02. Data Cleaning & Preparation:- 


*checking columns of the dataset:-* 


bombari_df.columns 


Index(['S#', 'Date', 
"City", 


"Killed Min', 


"Islamic Date', 
"Latitude', 
"Location Category', 
"Influencing Event/Event', 
"Killed Max", 
"No. of Suicide Blasts', 'Explosive Weight (max)', 


‘Longitude’, 'Province', 
"Location Sensitivity', 
"Target Type', 
"Injured Min', 


'Temperature(C)', 'Temperature(F)'], 


dtype='object') 


*checking basic information of the dataset:-* 


bombari_df.info() 


<class 'pandas.core.frame.DataFrame'> 
RangeIndex: 496 entries, Ø to 495 
Data columns (total 26 columns): 


dtypes: float64(7), int64(1), 
memory usage: 


f Column 


Date 

Islamic Date 
Blast Day Type 
Holiday Type 

Time 

City 

Latitude 
Longitude 
Province 

Location 

Location Category 
Location Sensitivity 
Open/Closed Space 


WONANEREWNRO 


BR 
RO 


Target Type 
Targeted Sect if any 
Killed Min 
Killed Max 
Injured Min 
Injured Max 
No. 


Hospital Names 
Temperature(C) 
25 Temperature(F) 


100.9+ KB 


Influencing Event/Event 


of Suicide Blasts 
Explosive Weight (max) 


Non-Null Count Dtype 
496 non-null int64 
496 non-null object 
342 non-null object 
486 non-null object 
72 non-null object 
285 non-null object 
496 non-null object 
493 non-null float64 
493 non-null object 
496 non-null object 
493 non-null object 
461 non-null object 
468 non-null object 
461 non-null object 
191 non-null object 
478 non-null object 
97 non-null object 
358 non-null float64 
480 non-null float64 
365 non-null float64 
464 non-null object 
414 non-null float64 
172 non-null object 
297 non-null object 
491 non-null float64 
489 non-null float64 
object(18) 


*correcting the data types of the columns:-* 


bombari df['Killed Min'] 
bombari df['Killed Max'] 


"Blast Day Type", 


"Location", 


"Holiday Type", 


Latitude 


33.718000 


24.991800 


24.991800 


24.991800 


30.209500 


28.571051 


NaN 


30.221057 


33.970623 


30.139626 


"Time", 


"Open/Closed Space", 


bombari df['Killed Min'].astype('Int64') 
bombari df['Killed Max'].astype('Int64') 


'Targeted Sect if any', 
'Injured Max', 
'Hospital Names', 


Longitude 


73.0718 


66.9911 


66.9911 


66.9911 


67.0182 


67.496895 


NaN 


67.002524 


71.43862 


66.981476 


Province 


Capital 


Sindh 


Sindh 


Sindh 


Baluchistan 


Balochistan 


Balochistan 


Balochistan 


KPK 


Balochistan 


Targeted 
Sect if 
any 


NaN 


NaN 


Christian 


Christian 


Shiite 


NaN 


NaN 


NaN 


NaN 


NaN 


Killed 
Min 


14.0 


NaN 


13.0 


NaN 


21.0 


7.0 


2.0 


1.0 


4.0 


Killed 
Max 


3.0 


15.0 


47.0 


22.0 


TO 


2.0 


3.0 


4.0 


Injured 
Min 


NaN 


NaN 


20.0 


NaN 


NaN 


22.0 


8.0 


6.0 


16.0 


Injured 
Max 


60 


40 


51 


65 


NaN 


22 


22 


No. of Explosive 


Suicide 
Blasts 


2.0 


1.0 


1.0 


1.0 


1.0 


1.0 


1.0 


1.0 


1.0 


1.0 


Weight 
(max) 


NaN 


NaN 


2.5 Kg 


NaN 


NaN 


NaN 


NaN 


15KG 


20KG 


10KG 


Hospital Names 


NaN 


NaN 


1Jinnah 
Postgraduate 
Medical Center 2. 
Civil ... 


NaN 


1.CMH Quetta 
\n2.Civil Hospital 3. 
Boland Medi... 


Larkana Hospital, 
Gandawah District 
Headquarte... 


NaN 


Combined Military 
Hospital (CMH)- 
(Quetta), Civ... 


Hayatabad Medical 
Complex(Peshawar) 


Civil 
Hospital(Quetta) 


bombari_df['Injured Min'] = bombari df['Injured Min'].astype('Int64') 

bombari df['Injured Max'] = pd.to numeric(bombari df["Injured Max"], errors= 'coerce') 
bombari df['Injured Max'] = bombari df['Injured Max'].astype('Int64') 

bombari df['No. of Suicide Blasts'] - bombari df['No. of Suicide Blasts'].astype('Int64') 


*getting only 'years' from the ‘date column’ Because they can be more useful in visualizations:-* 


bombari df['Year'] = bombari df.Date.str[-4:] 


*Extracting ‘Islamic Months' from the ‘Islamic Date’ column Because they can be more useful in visualizations:-* 


bombari df['Islamic Date'].isnull().sum() 


154 
bombari df['Islamic Date'] = bombari df['Islamic Date'].fillna('') 


bombari df["Islamic Date"] - bombari df["Islamic Date"].replace("","date unknown data month islamic year A.H") 
bombari df["Islamic Months"] - bombari df["Islamic Date"].apply(lambda x : x.split()[1]) 


bombari df['Islamic Months'].unique() 


array(['Jumaada', 'SHa'baan', 'safar', 'Raby'', 'Thw', 'MuHarram', 
'unknown data', 'Rajab', 'SHawwal', 'RamaDHaan', 'Shaban', 'Zeqad', 
'Zilhaj', 'Jamadi-ul-Awal', 'Jumada', 'Shawaal', 'Safar', 'Rabi'], 
dtype=object) 


*replacing with the correct islamic month names:-* 


bombari df["Islamic Months"] = bombari_df["Islamic Months"].replace( ("Jumaada" : "Jumaada al-Thaany", "Jumada": "Jumaada al-Thaany" 
"MuHarram": "Muharram", "safar": "Safar", "Rabi": "Rabi-ul-Awwal", "Raby": "Rabi-ul-Thaany", 
"Jamadi-ul-Awal": "Jamadi-ul-Awwal", "SHabaan": "Shaban", 'SHa'baan': 'Shaban', "RamaDHaan" : "Ramazan", 
"SHawwal": "Swawwal", "Thw": "Zilhaj")) 


bombari df['Islamic Months'].unique() 


array(['Jumaada al-Thaany', 'Shaban', 'Safar', 'Raby^', 'Zilhaj', 
'Muharram', 'unknown data', 'Rajab', 'Swawwal', 'Ramazan', 'Zeqad', 
'Jamadi-ul-Awwal', 'Shawaal', 'Rabi-ul-Awwal'], dtype=object) 


*checking 'Blast Day Type' column:-* 


bombari df['Blast Day Type'].unique() 
array(['Holiday', 'Working Day', nan, 'Weekend'], dtype-object) 


*checking 'Holiday Type' column:-* 


bombari df['Holiday Type'].unique() 


array(['Weekend', nan, 'Christmas/birthday of Quaid-e-Azam', 'Ashura', 
'Eid Milad un-Nabi', 'Iqbal Day', 'Eid-ul-azha', 'Labour Day', 
'Eid-ul-Fitar', 'Pakistan Day', 'Defence Day', 
'Christmas/ birthday of Quaid-e-Azam', 'General Elections' 
'Eid Holidays', 'Ashura Holiday', 'Eid ul Azha Holiday'], 
dtype=object) 


*making data more accurate:-* 


bombari df['Holiday Type'] = bombari df['Holiday Type'].replace(('Christmas/ birthday of Quaid-e-Azam': 'Christmas/birthday of Quaid-e-Azam', 
'Ashura': ‘Ashura Holiday', 'Eid ul Azha Holiday': 'Eid-ul-azha', 'Eid-ul-azha': 'Eid Holidays', 'Eid-ul-Fitar':'Eid Holidays", 
"Eid-ul-azha': 'Eid Holidays')) 


bombari df['Holiday Type'].unigue() 


array(['Weekend', nan, 'Christmas/birthday of Quaid-e-Azam', 
'Ashura Holiday', 'Eid Milad un-Nabi', 'Iqbal Day', 'Eid Holidays', 
'Labour Day', 'Pakistan Day', 'Defence Day', 'General Elections' 
'Eid-ul-azha'], dtype=object) 


*dealing with the 'City' Column:-* 


bombari df['City'].unique() 


array(['Islamabad', 'Karachi', 'Karachi ', 'Quetta', 'Rawalpindi', 
'North waziristan', 'Kohat', 'Attock ', 'Sialkot', 'Lahore ', 
'Swat', 'Hangu', 'Bannu', 'North Waziristan', 'Lasbela', 
'Malakand ', 'Peshawar', 'Peshawar ', 'D.I Khan ', 'Lakki Marwat', 
'Tank', ‘Islamabad ', 'Quetta ', 'Gujrat', 'Charsadda', 'Tank ' 
'Swat ', 'D.I Khan', 'Kohat ', 'Kuram Agency ', 'Shangla ' 
'Bajaur Agency', 'South Waziristan', 'Haripur', 'bannu', 
'karachi ', 'Sargodha', 'Nowshehra', 'Charsadda ', 
'Mohmand agency', 'Dara Adam Khel', 'Khyber Agency', 'Mardan', 
'Bhakkar ', 'Orakzai Agency', 'Mohmand Agency', 'Buner', 
'D.G Khan ', 'Pishin', 'Chakwal ', 'Lahore', 'Upper Dir', 
'Muzaffarabad', 'swat', 'Totalai', 'Multan', 'D.G Khan', 
'Lower Dir', 'Lakki marwat', 'Sudhanoti', 'Khyber Agency ' 
'peshawar', 'Poonch', 'Mansehra', 'Karak', 'Mohmand Agency ', 
'South waziristan', 'lakki marwat', 'Swabi ', 'Dara Adam khel', 
'karachi', 'Shikarpur', 'Sukkur', 'Chaman', 'D. I Khan', 
'Khanewal', 'Fateh Jang', 'Taftan', 'Tirah Valley', 'Wagah', 
'Zhob', 'KURRAM AGENCY', 'ATTOCK', 'MULTAN', 'Taunsa', 'Jacobabad', 
'Shabqadar-Charsadda', 'Khuzdar', 'Ghallanai, Mohmand Agency ' 
'Hayatabad', 'Mosal Kor, Mohmand Agency', 'Sehwan town', 


"Tangi, Charsadda District'], dtype=object) 


*making data consistent & accurate:-* 


bombari df['City'] = bombari df['City'].str.strip().str.lower() 
cities - bombari df['City'].unique() 
# cities 


4 say no to warnings... 
import warnings 
warnings.filterwarnings('ignore') 


import numpy as np 
from fuzzywuzzy import fuzz 


# Create an empty List to store the similar strings 
similar strings - [] 


# Loop through each city in the array 
for city in cities: 


# Loop through the rest of the cities 
for other city in cities: 
# Calculate the Levenshtein distance between the two strings 
distance = fuzz.ratio(city, other_city) 
# If the distance is greater than or equal to 8@ (i.e. they are very similar), and they are not the same string 
if distance >= 8@ and city != other_city: 
# Add the pair of similar strings to the List 
if sorted((city, other_city)) not in similar_strings: 
similar_strings.append(sorted((city, other_city))) 


similar_strings 


[[ "north waziristan', ‘south waziristan'], 
['d.g khan', 'd.i khan'], 

['d. i khan', 'd.i khan'], 

['kuram agency', 'kurram agency'], 

['d. i khan', 'd.g khan']] 


*replacing with the correct City names:-* 


bombari_df.City = bombari df.City.replace(('d. i khan': 'd.i khan", 'kuram agency': 'kurram agency')) 


*Results:-* 


bombari df.City.unique() 


array(['islamabad', 'karachi', 'guetta', 'rawalpindi', ‘north waziristan', 
'kohat', 'attock', 'sialkot', 'lahore', 'swat', 'hangu', 'bannu', 
'lasbela', 'malakand', 'peshawar', 'd.i khan', 'lakki marwat', 
'tank', 'gujrat', 'charsadda', 'kurram agency', 'shangla', 
'bajaur agency', 'south waziristan', 'haripur', 'sargodha', 
'nowshehra', 'mohmand agency', 'dara adam khel', 'khyber agency', 
'mardan', 'bhakkar', 'orakzai agency', 'buner', 'd.g khan', 
'pishin', 'chakwal', 'upper dir', 'muzaffarabad', 'totalai', 
'multan', 'lower dir', 'sudhanoti', 'poonch', 'mansehra', 'karak', 
'swabi', 'shikarpur', 'sukkur', 'chaman', 'khanewal', 'fateh jang', 
'taftan', 'tirah valley', 'wagah', 'zhob', 'taunsa', 'jacobabad', 
'shabqadar-charsadda', 'khuzdar', 'ghallanai, mohmand agency', 
'hayatabad', 'mosal kor, mohmand agency', 'sehwan town', 
'tangi, charsadda district'], dtype=object) 


*checking 'Province' column:-* 


bombari df.Province.unique() 


array(['Capital', 'Sindh', 'Baluchistan', 'Punjab', 'FATA', 'KPK', 'AJK', 
'Fata', 'Balochistan'], dtype-object) 


*replacing with the correct names:-* 


bombari df.Province = bombari df.Province.replace(('Baluchistan': 'Balochistan', 'Fata': 'FATA'}) 
bombari df.Province.unique() 


array(['Capital', 'Sindh', 'Balochistan', 'Punjab', 'FATA', 'KPK', 'AJK'], 
dtype=object) 


*checking 'Location' column:-* 


bombari_df.Location.sample(15) 


436 Imambargah 
398 Bolan Medical Complex 
298 WFP office Khar 
405 Torkhan customs office in Landikotal 
48 Training Center in Hangu 
298 Batakhel area 
368 Two suicide bombers blew themselves up at a cr... 
269 FC chowk near Dean's trade center 
146 near a check-post in Officers’ Colony 
160 Mosque in Jamrood Tehsil on peshawar-turkham h... 
16 Shrine of Bari Imam 
284 Mohmad Agency-Ghalani area-Political \nAgent 0... 
305 Near DCO office-Nowshehra road 
181 Quetta-Khi RCD highway-town of Torkham 
262 Bilalabad area near Khar town 


Name: Location, dtype: object 


*checking 'Location Category' column:-* 


bombari_df['Location Category'].unigue() 


array(['Foreign', 'Office Building', 'Hotel', 'Religious', 'Mobile', 
'Military', 'Residence', 'Park/Ground', 'Market', 'Police', 
'Airport', 'Government', 'Transport', 'Hospital', 'Educational', 
'Civilian', 'Bank', 'Foreigner', 'Government/Office Building', 
'Government Official', 'foreign', nan, 'Residential Building', 
'Commercial/residence', ' ', 'Highway'], dtype-object) 


*making data accurate:-* 


bombari df['Location Category'] = bombari df['Location Category'].replace(('Foreigner': 'Foreign', 'foreign': 'Foreign', 
'Government/Office Building': 'Office Building', 'Residential Building': 'Residence', ' ': np.nanj) 


bombari df['Location Category'].unique() 


array(['Foreign', 'Office Building', 'Hotel', 'Religious', 'Mobile', 
'Military', 'Residence', 'Park/Ground', 'Market', 'Police', 
'Airport', 'Government', 'Transport', 'Hospital', 'Educational', 
'Civilian', 'Bank', 'Government Official', nan, 
'Commercial/residence', 'Highway'], dtype-object) 


*overview of all columns:-* 


bombari df.columns 


Index(['S#', 'Date', 'Islamic Date', 'Blast Day Type', 'Holiday Type', 'Time', 
'City', 'Latitude', 'Longitude', 'Province', 'Location', 
'Location Category', 'Location Sensitivity', 'Open/Closed Space', 
'Influencing Event/Event', 'Target Type', 'Targeted Sect if any', 
'Killed Min', 'Killed Max', 'Injured Min', 'Injured Max', 
'No. of Suicide Blasts', 'Explosive Weight (max)', 'Hospital Names', 
'Temperature(C)', 'Temperature(F)', 'Year', 'Islamic Months'], 
dtype='object') 


*checking ‘Location Sensitivity’ column:-* 


bombari df['Location Sensitivity’ ].unigue() 
array(['High', 'Low', 'Medium', nan, 'low'], dtype=object) 


*making data accurate:-* 


bombari df['Location Sensitivity'] = bombari df['Location Sensitivity'].replace(('low': 'Low')) 
bombari df['Location Sensitivity'].unique() 
array(['High', 'Low', 'Medium', nan], dtype=object) 


*checking 'Open/Closed Space' column:-* 


bombari df['Open/Closed Space'].unique() 


array(['Closed', 'Open', 'open', 'Open 
dtype=object) 


» nan, 'closed', 'Open/Closed'], 


*making data accurate:-* 


bombari df['Open/Closed Space'] = bombari df[ 'Open/Closed Space'].str.strip() 
bombari df['Open/Closed Space'] = bombari df['Open/Closed Space'].replace(('open': 'Open', 'closed': 


bombari df['Open/Closed Space'].unique() 
array(['Closed', 'Open', nan, 'Open/Closed'], dtype-object) 


*checking 'Influencing Event/Event' column:-* 


bombari df['Influencing Event/Event'].sample(15) 


301 NaN 
473 Attack on Dargah Shah Noorani in Khuzdar district 
166 NaN 
253 Leader of anti-militants lashkar 
227 NaN 
249 NaN 
356 NaN 
459 Attack on Polio Centre 
404 It is said that Bannu TTP Militant commander S... 
484 Terrorist ambushed a police party while it was... 
478 inauguration of the Out Patient Department (OP... 
254 NaN 
313 NaN 
465 NaN 
12 Friday prayer 


Name: Influencing Event/Event, dtype: object 


*checking 'Target Type' column:-* 


bombari df['Target Type'].unigue() 


array(['Foreigner', 'Media', 'Religious', 'Military', 
'Government Official', 'civilian', 'Civilian', 'Police', 
'Government official', 'police', 'Children/Women', 
'Anti-Militants', 'foreigner', 'Unknown', 'religious', nan, 
'Rangers', 'Police & Rangers', 'Civilian & Police', 'Army' 
'Frontier Corps ', 'advocates (lawyers)', 'Civilian Judges', 
'Shia sect', 'Judges & lawyers'], dtype-object) 


*making data more accurate:-* 


bombari df['Target Type'] = bombari df['Target Type'].str.lower().str.strip() 
bombari df['Target Type'] - bombari df['Target Type'].replace(('civilian judges': 'judges & lawyers', 
'police': 'police & rangers', 'rangers': 'police & rangers')) 


bombari df['Target Type'].unique() 


array(['foreigner', 'media', 'religious', 'military', 
'government official', 'civilian', 'police & rangers', 
'children/women', 'anti-militants', 'unknown', nan 
'civilian & police', 'army', 'frontier corps', 'judges & lawyers', 
'shia sect'], dtype-object) 


*checking ‘Targeted Sect if any’ column:-* 


bombari df['Targeted Sect if any'].unique() 


array([nan, 'Christian', 'Shiite', 'shiite', 'Shiite/sunni', 'Sunni', 
"Jews', 'Ahmedi'], dtype=object) 


*making data more accurate:-* 


bombari df['Targeted Sect if any'] = bombari df['Targeted Sect if any'].replace(('shiite': 'Shiite')) 
bombari df['Targeted Sect if any'].unique() 


array([nan, 'Christian', 'Shiite', 'Shiite/sunni', 'Sunni', 'Jews', 
'Ahmedi'], dtype=object) 


*checking 'Explosive Weight (max)' column:-* 


bombari df['Explosive Weight (max)'].unique() 


array([nan, '2.5 Kg', '3@kg in each car', '2kg', '10-15 kg', '5 kg', 
'10 kg ', '5 to 6 Kg', '10 to 15 kg', 'more than 5Kg', '5 to 6 kg" 
'15 kg', '3 kg explosive and 3kg ball bearings', '6-7kg', '4-5kg' 
'SOkg', '30kg ', '100Kg', '150Kg', '1000 Kg', ' 8 Kg', '6G0Kg', 
'8 ft deep crater', '40 kg', '7 or 8Kg', ' 10 kg ', '8Kg', '100kg" 
'100', '85', '125 kg', '5kg', '160', '500 Kg', '24Kg', '10Kg', 
'20 Kg', '12Kg', '12 Kg', '150 kg ', '180Kg', '100 Kg', '70 kg ', 
'16 Kg', '15Kg', '10 Kg', '7Kg ', '400Kg', '200Kg', '50-60Kg', 
'250Kg', '1000 KG', '250 kg', '600 kg', '40Kg', '5Kg', ' 12 kg', 
'600Kg', '12 kg', '14 kg', '8 kg', '300 kg', '20 kg ', '10 kg', 
'20Kg', ' 6 kg', '1200Kg', '12 kg ', '1000 Kg ', '400 KG', 
'1000 kg', '7Kg', '9Kg', '16Kg', '250kg', 'Skg', '180kg', 
'16-20kg', '300kg', '450kg', '7kg', '6Økg', '500 kg', '5-6kg', 
'3kg', '6kg', '6Kg', '100-110kg', '30kg', '7kg-8kg', '100kg-150kg', 
'4kg-5kg', '12kg', '600kg-800kg', '4kg-6kg', '1000kg', '7-8kg', 
'6-7Kg', '20-25kg', '25-30kg', '14kg', '10kg', '45KG', '15KG', 
'13KG', '8KG', '8-10KG', '6KG', '8-10 KG', ' 20 KG', '15 KG', 
'75KG', '25KG', '20KG', '10KG'], dtype-object) 


*making data consistent & accurate:-* 


"Closed')) 


"advocates (lawyers)': 


‘judges € lawyers’, 


bombari df['Explosive Weight (max)'] = bombari df[ 'Explosive Weight (max)'].str.strip() 

bombari df['Explosive Weight (max)'] = bombari df[ 'Explosive Weight (max)'].str.replace('KG', 'kg') 
bombari df['Explosive Weight (max)'] = bombari df[ 'Explosive Weight (max)'].str.replace('Kg', 'kg') 
bombari df['Explosive Weight (max)'] = bombari df[ 'Explosive Weight (max)'].str.replace(' ', '') 


bombari df['Explosive Weight (max)'] = bombari df[ 'Explosive Weight (max)'].replace({'3@kgineachcar': '30kg', '3kgexplosiveand3kgballbearings': 
'morethan5kg 'ekg', '7kg-8kg': '7-8kg', '100kg-150kg': '100-150kg', ‘4kg-5kg': '4-5kg', '600kg-800kg': '600-800kg', 
'Akg-6kg': '4-6kg', '5to6kg': '6kg', '10tol5kg': '12kg', '7or8kg': '7kg')) 


bombari df['Explosive 


Weight (max)'].replace('8ftdeepcrater', np.nan, inplace-True) 


# bombari df['Explosive Weight (max)'].unique() 


result - 
result 


['10-15kg', 
'6-7kg', 
'4-5kg', 
'50-60kg' , 
'16-20kg', 
'5-6kg', 
'100-110kg', 
"7-8kg', 
'100-150kg', 
'600-800kg' , 
'4-6kg', 
'20-25kg', 
"25-30kg', 
"8-1ekg' ] 


bombari df['Explosive Weight (max)'] = bombari df[ 'Explosive Weight 


'10-15kg': '12kg', 
'6-7kg': '6kg', 
'4-5kg': 'Akg', 
'50-60kg': '55kg', 
'16-20kg': '18kg', 
'5-6kg': '5kg', 
'100-110kg': '105kg', 
'7-8kg': '7kg', 
'100-150kg': '125kg', 
'600-800kg': '700kg', 
'A-6kg': '5kg', 
'20-25kg': '22kg', 
'25-30kg': '27kg', 
'8-10kg': '9kg" 

}) 

*Results:-* 


bombari_df['Explosive Weight (max)'].unique() 


[s for s in bombari df['Explosive Weight (max)'].unique() if '-' in str(s)] 


(max) ' ]. replace(( 


array([nan, '2.5kg', '30kg', '2kg', '12kg', '5kg', '10kg', '6kg', '15kg', 
'Akg', 'Søkg', '100kg', '150kg', '1000kg', '8kg', '60kg', '4Əkg', 
'7kg', '100', '85', '125kg', '160', '500kg', '24kg', '20kg' 
'180kg', '70kg', '16kg', '400kg', '200kg', '55kg', '250kg', 
'600kg', '14kg', '300kg', '1200kg', '9kg', '18kg', '450kg', '3kg' 
'105kg', '700kg', '22kg', '27kg', '45kg', '13kg', '75kg', '25kg'] 


dtype=object) 


*checking 'Temperature(C)' colummn:-* 


bombari df['Temperature(C)'] 


e 15.835 
1 23.770 
2 31.460 
3 31.430 
4 33.120 
491 27.800 
492 26.000 
493 25.000 
494 22.000 
495 29.000 


Name: Temperature(C), Length: 496, dtype: float64 


*making data accurate:-* 


bombari df = bombari df.replace([np.inf, -np.inf], np.nan).dropna(subset-['Temperature(C)']) 


bombari df['Temperature(C)'] - bombari df['Temperature(C)'].round(). 


bombari df['Temperature(C)'] 


e 16 
1 24 
2 31 
3 31 
4 33 
491 28 
492 26 
493 25 
494 22 
495 29 


Name: Temperature(C), Length: 491, dtype: int32 


*making 'Hospital Names' column data more consistent:-* 


bombari df['Hospital Names'].head(10) 


NaN 
NaN 
1.Jinnah Postgraduate Medical Center 2. Civil ... 
NaN 
1.CMH Quetta \n2.Civil Hospital 3. Boland Medi... 
1.District headquarters \nHospital 
NaN 
Civil hospital 
Civil hospital-Liaquat National hospital-Jinna... 
NaN 


YVAØNAUSRUNKO 


Name: Hospital Names, dtype: object 


*let's look at our dataset now:-* 


bombari df 


astype(int) 


'ékg', 


491 


492 


493 


494 


495 


Islamic 

SH Date Date 
Sunday- 22 
Jumaada 

1 November al-THaany 
19195 1416 AH 
Monday- 10 


2 November  SHa'baan 
6-2000 1421 AH 


Wednesday- 25 safar 
May 8-2002 1423 AH 


Friday-June lis) dle 
14-2002  THaany 
1423 A.H 

Friday-July ae 
al-awal 

49035 4424 AH 
Thursday- 13 


492 October-05- MuHarram 
2017 1439 AH 


Thursday- 27 


493 October-19- MuHarram 

2017 1439 AH 

Thursday- 19 Safar 

494  November- 1439 AH 
9-2017 i 


Friday- — 5 Rabi Al- 
495 November- Awwal 
24-2017 1439 AH 


Saturday- 6 Rabi Al- 
496 November- Awwal 
25-2017 1439 A.H 


491 rows x 28 columns 


Blast 
Day 
Type 


Holiday 


Working 
Day 


Working 
Day 


Working 
Day 


Working 


Day 


Working 


Day 


Working 
Day 


Working 
Day 


Working 
Day 


Working 
Day 


Holiday 
Type 


Weekend 


NaN 


NaN 


NaN 


NaN 


NaN 


NaN 


NaN 


NaN 


NaN 


Time City 


NaN islamabad 


NaN karachi 


7:45 AM karachi 


11:10:00 


AM karachi 


NaN quetta 


NaN quetta 


8:00 AM quetta 


7:00 AM quetta 


7:00AM peshawar 


7:00 AM quetta 


03. Exploratory Analysis & Visualizations:- 


import matplotlib.pyplot as plt 


import seaborn as sns 
import plotly.express as px 


*columns overview:-* 


bombari_df.columns 


Index(['S#', 'Date', 'Islamic Date', 


'City', 'Latitude', 'Longitude', 


dtype='object') 


"Blast Day Type', 
"Province", ‘Location’, 
"Location Category’, ‘Location Sensitivity’, 
"Influencing Event/Event', 
"Killed Min', ‘Killed Max", 
"No. of Suicide Blasts', ‘Explosive Weight (max)', 
"Temperature(C)', 'Temperature(F)', 'Year', 


"Target Type', 
"Injured Min', 


01. Top 10 Cities with Most Blasts:- 


*plotting:-* 


‘Injured Max', 


"Holiday Type', 


Latitude 


33.718000 


24.991800 


24.991800 


24.991800 


30.209500 


28.571051 


NaN 


30.221057 


33.970623 


30.139626 


"Open/Closed Space', 
‘Targeted Sect if any", 


"Hospital Names", 


"Islamic Months'], 


plot 1 = bombari df['City'].value counts().head(10) 


pit. 


figure(figsize-(16, 7)) 


plot 1.plot(kind-'bar', color-'teal') 


plt 


for 


plt. 


pit. 


xlabel('Cities', fontsize-19) 


ylabel('No. of Blasts', fontsize-18) 


i in range(len(plot 1)): 


plt.text(i, plot 1[i]+1, plot 1[i], ha-'center', fontsize=14) 


.title('Top 10 Cities with Most Blasts', fontsize-20, fontweight-'bold') 
plt. 
pits 


xticks(range(len(plot 1)), plot 1.index, fontsize=14, fontweight-'bold') 


show() 


'Time', 


Longitude 


73.0718 


66.9911 


66.9911 


66.9911 


67.0182 


67.496895 


NaN 


67.002524 


71.43862 


66.981476 


Province 


Capital 


Sindh 


Sindh 


Sindh 


Balochistan 


Balochistan 


Balochistan 


Balochistan 


KPK 


Balochistan 


Killed 
Max 


47 


22 


Injured 
Min 


«NA» 


«NA» 


20 


«NA» 


«NA» 


22 


16 


Injured 


Nex Sie 
60 2 

3 1 

40 1 

51 1 

65 1 
<NA> 1 
22 1 

8 1 

8 1 

22 1 


No. of Explosive 


Weight 
(max) 


NaN 


NaN 


2.5kg 


NaN 


NaN 


NaN 


NaN 


15kg 


20kg 


10kg 


Hospital Names 


NaN 


NaN 


1Jinnah 
Postgraduate 
Medical Center 2. 
Civil ... 


NaN 


1.CMH Quetta 
\n2.Civil Hospital 3. 
Boland Medi... 


Larkana Hospital, 
Gandawah District 
Headquarte... 


NaN 


Combined Military 
Hospital (CMH)- 
(Quetta), Civ... 


Hayatabad Medical 
Complex(Peshawar) 


Civil 
Hospital(Quetta) 


Temperature(C) T 


16 


24 


31 


31 


33 


28 


26 


25 


22 


29 


In [ 


In [ ] 


Top 10 Cities with Most Blasts 


707 
604 
50 4 


404 


75 
37 
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204 
104 
0 
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02. Sucide Blasts Locations in Pakistan:- 


*extracting pakistan's latitudes 8. longitudes for the map:-* 

import geopandas as gpd 

from shapely.geometry import mapping 

world = gpd.read file(gpd.datasets.get path( 'naturalearth lowres')) 


poly = world[world.name == "Pakistan"].geometry 
poly_mapped = mapping(poly) 


poly coordinates = poly mapped["features"][0]['geometry' ][ 'coordinates'][0] 
poly dict = [('lat': coords[1],'lon': coords[0]) for coords in poly coordinates] 
pak boundary coord - pd.DataFrame(poly dict) 


pak boundary coord 


lat lon 
0 35494010 77.837451 
1 34.653544 76.871722 
2 34.504923 75.757061 
3 34.748887 74.240203 


4 34317699 73.749948 


61 37.020841 74.575893 
62 37.133031 75.158028 
63 36.666806 75.896897 
64 35.898403 76.192848 


65 35494010 77.837451 


66 rows x 2 columns 
*overview:-* 


# fig = px.line(pak boundary coord, x="Lon", y-"Lat", height-800, width-1300) 
4 fig.update_traces(Line_color="darkgreen") 


*checking city names:-* 


bombari df.City.unique() 


array(['islamabad', 'karachi', 'quetta', 'rawalpindi', 'north waziristan', 
'kohat', 'attock', 'sialkot', 'lahore', 'swat', 'hangu', 'bannu', 
'lasbela', 'malakand', 'peshawar', 'd.i khan', 'lakki marwat', 
'tank', 'gujrat', 'charsadda', 'kurram agency', 'shangla', 
'bajaur agency', 'south waziristan', 'haripur', 'sargodha', 
'nowshehra', 'mohmand agency', 'dara adam khel', 'khyber agency', 
'mardan', 'bhakkar', 'orakzai agency', 'buner', 'd.g khan', 
'pishin', 'chakwal', 'upper dir', 'muzaffarabad', 'totalai', 
'multan', 'lower dir', 'sudhanoti', 'poonch', 'mansehra', 'karak', 
'swabi', 'shikarpur', 'sukkur', 'chaman', 'khanewal', 'fateh jang', 
'taftan', 'tirah valley', 'wagah', 'zhob', 'taunsa', 'jacobabad', 
'shabqadar-charsadda', 'ghallanai, mohmand agency', 'hayatabad', 
'sehwan town', 'tangi, charsadda district'], dtype=object) 


*feeding correct latitudes and longitudes for the above cities in the dataset:-* 


city dict - ('islamabad': (33.6844, 73.0479), 
'karachi': (24.8607, 67.0011), 


north wazi 


Cities 


'quetta': (30.1798, 66.9750), 
'rawalpindi': (33.5651, 73.0169), 
'north waziristan': (32.9881, 69.9344), 
'"kohat': (33.5869, 71.4417), 
Kattock4 (33.7717, 72.3592), 
'sialkot': (32.4947, 74.5319), 
'lahore': (31.5820, 74.3294), 

'swat': (34.7566, 72.3558), 

'hangu': (33.5292, 71.0728), 

'bannu': (32.9869, 70.6050), 
'lasbela': (26.8531, 66.3833), 
'malakand': (34.5500, 71.9167), 
'peshawar': (34.0150, 71.5250), 

'd.i khan': (31.8297, 70.9017), 
'lakki marwat': (32.6144, 70.9125), 
'tank': (32.2219, 70.3797), 

'gujrat': (32.5722, 74.0750), 
'charsadda': (34.1450, 71.7406) 
'kurram agency': (33.8813, 69.8922), 
'shangla': (34.7953, 72.4753), 
'bajaur agency': (34.8103, 71.8158), 
'south waziristan': (32.3617, 69.8553), 
'haripur': (33.9990, 72.9343), 
'sargodha': (32.0836, 72.6711), 
'nowshehra': (34.0158, 71.9819), 
'mohmand agency': (34.1658, 71.9244), 
'dara adam khel': (33.8439, 71.8864), 
'khyber agency': (34.0511, 71.1083), 
"mardan": (34.1987, 72.0499), 
"bhakkar": (31.6267, 71.0656), 
"orakzai agency": (33.1025, 70.1394), 
"buner": (34.5278, 72.4786), 

"d.g khan": (30.0561, 70.6342), 
"pishin": (30.5806, 66.9961), 
"chakwal": (32.9308, 72.8583), 
"upper dir": (35.2058, 71.8761), 
"muzaffarabad": (34.3700, 73.4711), 
"totalai": (34.6342, 72.6522), 
"multan": (30.1989, 71.4783), 
"lower dir": (34.8156, 71.8342), 
"sudhanoti": (33.7947, 73.5836), 
"poonch (33.7669, 74.1070), 
"mansehra": (34.3334, 73.1968), 
"karal (33.1167, 71.0833), 

"swabi": (34.1229, 72.4726), 
“shikarpur": (27.9555, 68.6375), 
"sukkur": (27.7136, 68.8486), 
(30.9186, 66.4464), 
(30.3012, 71.9329), 
"fateh jang": (33.3918, 72.9776), 
"taftan": (29.7979, 61:7027), 

"tirah valley": (33.7466, 70.1176), 
"wagah": (31.6030, 74.4403), 

"zhob": (31.3416, 69.4490), 

"taunsa": (30.6975, 70.6722), 
"jacobabad": (28.2810, 68.4390), 
"shabqadar-charsadda": (34.1354, 71.8096) 
"khuzdar": (27.7384, 66.6432), 
"ghallanai, mohmand agency": (34.2224, 71.5558), 
"hayatabad": (33.9902, 71.4623), 
"mosal kor, mohmand agency": (34.2639, 71.5397), 
"sehwan town": (26.4241, 67.8609), 
"tangi, charsadda district": (34.1835, 71.7525)) 


for index, row in bombari df.iterrows(): 
# get the name of the city for the current row 
city - row['City'] 
4 get the Latitude and Longitude for the current city from the city dict 
lat, long - city dict[city] 
# replace the Latitude and Longitude values in the current row with the Lat and Long values 
bombari df.at[index, 'Latitude'] - lat 
bombari df.at[index, 'Longitude'] - long 


*plotting:-* 


[ import plotly.graph objects as go 
import plotly.offline as pyo 


fig = px.scatter(bombari df, x-'Longitude', y='Latitude', hover name='City', hover data-['Date', 'Province']) 
fig.update traces(marker-dict(color-'red', size-6), hoverlabel-dict(font-dict(family-'Arial', size-18, color-'white'))) 
fig.update layout(mapbox style-"open-street-map") 

fig.update layout(mapbox zoom-5, mapbox center - ("lat": 30.3753, "lon": 69.3451)) 

fig.add trace(px.line(pak boundary coord, x-"lon", y-"lat").data[0]) 

fig.update traces(line color-'darkgreen') 

fig.update layout(height-1040, width-1400, plot bgcolor-'darkkhaki') 


pyo.plot(fig, filename-'my plot.html', auto open-True) 


fig.update traces(hovertemplate-"«b»Blast Location: %(hovertext)</b><br>" + 
"«b»Date:«/b» %{customdata[@]}<br>" + 
"«b»Province:«/b» X(customdata[1])«br»" + 
"<extra></extra>") 


fig.update layout( 


title={ 
"text "Sucide Blasts Locations in Pakistan', 
"font ('size': 30, 'color': 'black') 
Le 44-00 JA 
'xanchor': "center" 


+ 

xaxis title='Longitude', 

xaxis title font-('size': 22, 'color': 'black'], 
yaxis title font-('size': 22, 'color': 'black') 


display(fig) 


03. Average No. of People Killed & Injured per Year:- 


*dealing:-* 


bombari df['Average Killed'] - bombari df['Killed Max'] * bombari df['Killed Min'] / 2 
bombari df['Average Injured'] = bombari df['Injured Max'] + bombari df['Injured Min'] / 2 


bombari df['Average Killed'].fillna(0, inplace-True) 
bombari df['Average Injured'].fillna(0, inplace-True) 


70 


Longitude 


bombari df['Average Killed'] = bombari df['Average Killed'].astype('int') 
bombari df['Average Injured'] - bombari df['Average Injured'].astype('int') 


*grouping:-* 


result - bombari df.groupby("Year")[['Average Killed', 'Average Injured']].sum() 
result 


Average Killed Average Injured 


Year 

1995 22 0 
2000 0 0 
2002 21 50 
2003 95 0 
2004 102 458 
2005 105 284 
2006 227 129 
2007 978 2148 
2008 1146 2828 
2009 1224 4366 
2010 1500 3409 
2011 785 1757 
2012 301 875 
2013 917. 2258 
2014 267 643 
2015 265 360 
2016 453 1161 
2017 429 793 
*plotting:-* 


plt.figure(figsize-(14,6)) 


plt.plot(result.index, result['Average Killed'], label-'Killed', color-'red', linewidth-2) 
plt.plot(result.index, result['Average Injured'], label-'Injured', color-'blue', linewidth-2) 


plt.xlabel('Year', fontsize-17) 

plt.ylabel('Total No. of People', fontsize-15) 

plt.title('Average No. of People Killed & Injured per Year', fontsize-19, fontweight - 'bold') 
plt.legend(fontsize-12) 

plt.grid(alpha-0.3) 

plt.show() 
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04. Suicide Blasts w.r.t Islamic Months:- 


*plotting:-* 


import matplotlib.ticker as ticker 
plt.figure(figsize=(16, 6)) 


colors = sns.color palette('rocket') 
counts, edges, patches = plt.hist(bombari df['Islamic Months'], bins=12, rwidth=8.9, color=colors[1]) 


for count, edge, patch in zip(counts, edges, patches): 
if count > 0: 
plt.annotate(str(count), xy=(edge + 0.5, count), xytext-(0, 5), textcoords-'offset points', ha="center", va='bottom') 


plt.grid(axis-'x', alpha-0.75) 

plt.gca().set facecolor('4f2f2f2') 

plt.title('Suicide Blasts w.r.t Islamic Months', fontsize-20, fontweight-'bold', fontname-'Arial') 
plt.xlabel('Islamic Months', fontsize-17, fontname='Arial') 

plt.xticks(fontsize-14, fontweight-'bold', rotation-45, fontname-'Arial') 

plt.ylabel('No. of Suicide Blasts', fontsize-16, fontname-'Arial') 

plt.tick params(axis-'x', which-'major', pad-10) 
plt.gca().yaxis.set major formatter(ticker.FormatStrFormatter('%1.0f')) 

plt.subplots adjust(wspace-0.4) 

plt.legend(['Suicide Bombings'], loc-'upper left', fontsize-12) 


plt.show() 
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05. Percentage of Suicide Blasts By Location Sensitivity: 


*dealing:-* 


df loc sens = bombari_df.groupby([ ‘Location Sensitivity']).size().reset index(name-'Count') 


# Calculate the percentage of attacks for each Location sensitivity 
total attacks = df loc sens['Count'].sum() 
df loc sens['Percentage'] = df loc sens['Count'] / total attacks * 100 


*plotting:-* 


plt.figure(figsize- (17,6)) 


sns.set style("darkgrid") 
colors - ['orangered', 'lightsalmon', 'coral'] 


fig, ax = plt.subplots(figsize-(8, 6)) 
ax.pie(df loc sens['Percentage'], labels=df loc sens['Location Sensitivity'], autopct='%1.1f%%', colors-colors, wedgeprops-dict(width-0.6, edgecolor-'white'), textprops={'fontsi 


ax.set title('Percentage of Suicide Blasts by Location Sensitivity', fontsize-19, fontweight - 'bold') 
ax.legend(loc-'best', bbox to anchor-(1, 0.5)) 


circle = plt.Circle(xy-(0,0), radius-0.5, facecolor-'white') 
ax.add artist(circle) 


total attacks str = f"Total attacksWn(total attacks)" 
ax.text(0, 0, total attacks str, ha='center', va-'center', fontsize=16) 


ax.set xticks([]) 
ax.set yticks([]) 
ax.set xlabel('') 
ax.set ylabel('') 


plt.show() 


«Figure size 1700x600 with 0 Axes» 
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06. Total Victims w.r.t Target Type:- 


*dealing:-* 


bombari df['Total Victims'] = bombari df['Average Injured'] + bombari df['Average Killed'] 


*grouping:-* 
grouped data = bombari_df.groupby('Target Type')['Total Victims'].sum() 
*plotting:-* 


fig, ax = plt.subplots(figsize=(15, 6)) 
grouped data.plot(kind-'bar', stacked-True, ax=ax, color=colors) 


colors = ['orange', 'blue', 'gray', 'green', 'teal', '#FFA@7A', 'indianred', '#85144b',] 
sns.set style('whitegrid') 

ax.set xlabel('Target Type', fontsize-18) 

ax.set xticklabels(ax.get xticklabels(), fontdict-('fontsize': 12, 'fontweight': 'bold')) 
ax.set ylabel('Total Victims', fontsize-15) 

ax.set title('Total Victims w.r.t Target Type', fontsize-20, fontweight-'bold') 


for rect in ax.containers: 
ax.bar label(rect, label type-'edge', fontsize-11, fontweight='bold') 


Total Victims w.r.t Target Type 


ax.grid(axis-'y', linestyle="--") 
plt.show() 
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07. Suicide Blasts w.r.t Islamic Months:- 
*plotting:-* 


import random 


targeted sects = bombari df['Targeted Sect if any'].value counts() 
sizes - targeted sects.values.tolist() 


plt.figure(figsize-(10, 6)) 

plt.pie(sizes, colors-colors, startangle-90) 

plt.axis('equal') 

colors = ['£'«''.join(random.sample('0123456789abcdef', 6)) for i in range(len(sizes))] 


table colors - [plt.Rectangle((0, 0), 1, 1, fc-color) for color in colors] 


legend - plt.legend(table colors, targeted sects.index.tolist(), loc-'center left', bbox to anchor-(1, 0.5)) 


plt.setp(legend.get texts(), fontsize-'large') 
plt.title('Targeted Sectors in Suicide Blasts', fontsize-20, fontweight - 'bold') 


for label in legend.get texts(): 
label.set fontsize(14) 


print(targeted sects) 


plt.show() 
Targeted Sect if any 
Shiite 47 
Sunni 38 
Christian 9 
Shiite/sunni 1 
Jews 1 
Ahmedi 1 


Name: count, dtype: int64 
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08. No. of Suicide Blasts on Open 8 Closed Spaces:- 


*plotting*:-* 


space counts = bombari df.pivot table(index-'Year', columns='Open/Closed Space', values='S#', aggfunc-'count') 
plt.figure(figsize-(17,6)) 
sns.set style("whitegrid") 


ax = space counts.plot(kind-'bar', stacked-True, figsize-(12, 6), color=['41f77b4', 'iff7f0e'], edgecolor-'white', linewidth=0.5) 


ax.set title('No. of Suicide Blasts on Open 8 Closed Spaces', fontsize-20, fontweight='bold') 
ax.set xlabel('Year', fontsize-16) 

ax.set ylabel('Number of Attacks', fontsize-16) 

plt.xticks(fontsize-14, fontweight-'bold') 

plt.yticks(fontsize-14, fontweight-'bold') 


ax.legend(labels-['Closed Spaces', 'Open Spaces'], loc='upper left', fontsize-16, title-'Space Type') 
ax.grid(axis-'y', color='#666666', alpha=0.5) 


ax.spines['top'].set visible(False) 
ax.spines['right'].set visible(False) 


ax.axhline(0, color-'black', lw-1, alpha=0.5) 


for i in range(len(space counts)): 
closed count - space counts.iloc[i, 0] 
open count - space counts.iloc[i, 1] 
ax.text(i, closed count/2, format(closed count, '.0f'), ha='center', fontsize-14, fontweight-'bold', color='white') 
ax.text(i, closed count + open count/2, format(open count, '.0f'), ha='center', fontsize=14, fontweight="bold") 


plt.subplots adjust(hspace-0.4) 


plt.show() 


«Figure size 1700x600 with 0 Axes» 
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09. Temperature vs Suicide Blasts:- 


*plotting:-* 

plt.figure(figsize = (15,6)) 

plt.scatter(bombari df['Temperature(C)'], bombari df['St'], s=50, alpha-0.8, c-bombari df['Temperature(C)'], cmap='viridis') 
plt.xlabel('Temperature (C)', fontsize=15) 

plt.ylabel('No. of Suicide Blasts', fontsize=15) 


plt.title('Temperature vs Suicide Blasts', fontsize=20, fontweight = 'bold') 


cbar = plt.colorbar() 
cbar.set label('Temperature (C)', fontsize-15) 


plt.show() 


Temperature vs Suicide Blasts 


500 
e = eve 1 o op 


o 990 o ee 8 > ad ge 
o o. e eg e o 
*s 8,%,0 o 50 09 o 
400 . oe toos “of e. *$8 ES 
ee TLD 8 9 
oo e o e e 

o e e$ 8 vat eg “o 9° 00 9 
E e ° 4 0.2 o 
£n 300 e 
o 
3 
2 e 
5 
a 
« 200 
o 
o 
z 

100 e n 

e 
0 
0 10 20 30 40 


Temperature (C) 


8 


3 
Temperature (C) 


= 
> 


